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^ : Abstract 



A unique pattern selection in the absolutely unstable regime of a driven, 



> 

o 
o 

' nonlinear, open-flow system is analyzed: The spatiotemporal structures of 

' rotationally symmetric vortices that propagate downstream in the annulus 

^ ' of the rotating Taylor-Couette system due to an externally imposed axial 



through-flow are investigated for two different axial boundary conditions at 
the in- and outlet. Detailed quantitative results for the oscillation frequency, 



^ ■ the axial profile of the wave number, and the temporal Fourier amplitudes 

of the propagating vortex patterns obtained by numerical simulations of 
the Navier-Stokes equations are compared with results of the appropriate 
Ginzburg-Landau amplitude equation approximation and also with experi- 
ments. Unlike the stationary patterns in systems without through-fiow the 
spatiotemporal structures of propagating vortices are independent of parame- 
ter history, initial conditions, and system's length. They do, however, depend 
on the axial boundary conditions in addition to the driving rate of the inner 
cylinder and the through-flow rate. Our analysis of the amplitude equation 
shows that the pattern selection can be described by a nonlinear eigenvalue 
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problem with the frequency being the eigenvalue. The complex amplitude 
being the corresponding eigenfunction describes the axial structure of inten- 
sity and wave number. Small, but characteristic differences in the structural 
dynamics between the Navier-Stokes equations and the amplitude equation 
are mainly due to the different dispersion relations. Approaching the border 
between absolute and convective instability the eigenvalue problem becomes 
effectively linear and the selection mechanism approaches that one of linear 

front propagation. 
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I. INTRODUCTION 

In many nonlinear continuous systems dissipative structures branch out of a homoge- 
neous basic state when the external stress exceeds a critical threshold. Examples for these 
transitions are Taylor- Couette flow, Rayleigh-Benard convection, binary-fluid convection, 
flame-front propagation, and some chemical or biological processes Often, for a fixed 
configuration of parameters and boundary conditions a continuous or discrete family of pat- 
terns with different wave numbers is stable. Their stability regime, e. g., a band of wave 
numbers might be limited by the possibility of resonant triad interactions of modes like those 
described by the Eckhaus or Benjamin- Feir mechanism [|^. The stable structures within such 
a band can be generated by appropriately engineered time histories of the parameters and/or 
by properly changing the boundary conditions, e. g., the system size. The most intensively 
investigated examples in this respect are the structures of Taylor vortices [^-[^ in an annulus 
between concentric cylinders of which the inner one rotates and convective roll patterns in 
horizontal layers of one-component fluids [p],|8|-[Tl|,^ or binary mixtures heated from below 



12,13 



This multiplicity of solutions of the underlying nonlinear partial differential equations 
that stably coexist for a fixed configuration of parameters and boundary conditions seems 
to disappear in an open-fiow system: Recent numerical simulations of Rayleigh-Benard 
convective rolls traveling downstream in an imposed horizontal Poiseuille fiow showed [|lj-|l6 



that their structure is uniquely selected - i. e. it is independent of parameter history, initial 
conditions, and system size - in the absolutely unstable regime. This is the parameter 
regime of an open-fiow system in which the secondary pattern starting, e. g., from a spatially 
localized perturbation can grow in upstream as well as in downstream direction [O. By 



contrast, in the convectively unstable regime initial perturbations are blown out of 
the system - both, the upstream as well as the downstream facing front of the growing 
structure move downstream. In the absolutely unstable regime the structure expansion 
proceeds until the upstream (downstream) moving front encounters in a finite system the 
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inlet (outlet) and adjusts to the inlet (outlet) boundary condition. The final pattern resulting 
in such a situation shows a characteristic streamwise profile of the amplitude growing with 
increasing distance from the inlet and of the wave number variation, and a characteristic 
global oscillation frequency associated with the downstream motion of the pattern. 

In this work we elucidate in numerical and analytical detail how such a uniquely selected 
spatiotemporal pattern structure can be understood as a nonlinear eigenvalue problem with 
the oscillation frequency being the eigenvalue and the profiles of pattern intensity and wave 
number determining the corresponding eigenfunction. We also show how this pattern selec- 
tion process is related to the one occurring behind a front or "domain wall" that spatially 
separates an unstable, homogeneous state from a stable, structured state. To that end 
we present results of extensive numerical and analytical investigations of vortex patterns 
in the annulus of the Taylor-Couette setup with an externally imposed axial through-fiow. 
For small through-fiow rates and small rotation rates of the inner cylinder the structure 
of propagating vortex (PV) fiow with vortices being advected in downstream direction is 
rotationally symmetric. Only for higher through-fiow rates, which are not investigated here, 
there is a bifurcation [0,^] to spirals [pO|-p6| . We have performed numerical simulations of 



the full, rotationally symmetric, two dimensional (2D) Navier-Stokes equations (NSE). They 
are compared with our numerical and analytical results obtained from the appropriate ||2^ 



ID Ginzburg-Landau amplitude equation (GLE) approximation to the problem of PV fiow 
and with experimental resuhs [|g-|2|j2|j2|j2l 



The eigenvalue problem of the pattern selection by the imposed fiow can best be analyzed 
and explained within the GLE framework. Therein, the complex amplitude A{z, t) of the 
PV fiow depends in the absolutely unstable regime on the streamwise position z and on time 
t in a multiplicative way only, 

A{z,t) = a{z)e-'^\ (1.1) 

after transients have died out. The 2;-independent frequency fl and the complex, z-dependent 
amplitude a{z) that describes the streamwise variation of pattern modulus and wave number 
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are simultaneously fixed via a solvability condition in the form of a nonlinear eigenvalue 
problem. The selection of Q and a{z) seems to result from requiring the spatial variation of 
the amplitude to be as small as possible under the imposed boundary conditions for A at 
the ends of the annulus. 

Approaching the border between the absolutely and convectively unstable regimes the 
pattern selection mechanism becomes linear: For driving and through-flow rates on this 
border line the selected frequency is the one resulting from a linear front whose spatiotem- 
poral behavior is governed by the fastest growing linear mode. The latter is identified by a 
particular saddle of the complex linear dispersion relation over the complex wave number 
plane [Q. Now, the linear dispersion relations of NSE and GLE differ for supercritical 



control parameters. And therefore the PV structures selected by NSE or GLE differ in a 
characteristic way. 

It should be emphasized that the structural dynamics of pattern formation in the convec- 
tively unstable regime at larger through-flow rates and/or smaller driving rates substantially 
differs from the one investigated here in the absolutely unstable regime. The latter regime 
is governed by nonlinear contributions in the balance equations; the resulting patterns are 
uniquely selected and insensitive to initial conditions, parameter history, and small pertur- 
bations. On the other hand, in the convectively unstable regime that has attracted more 
experimental activities lately [p^ , ^ ,pl|-|35| the growing patterns are sensitive to initial con- 



ditions and perturbations. Thus, e. g., the noise sustained patterns [^] occurring in this 
regime depend on details of the spatiotemporal properties of the perturbation source. 



The Taylor- Couette system with an imposed axial through-flow has been investi- 
gated as a well deflned open-flow system theoretically and experimentally ||39| since 
the beginning 1930s. Linear stability analyses of the basic state to traveling axisymmetric 
vortices were performed using various approximations, e. g. , for the narrow-gap limit, or 
simplifled azimuthal velocity, or axial through-flow proflles EOl, and later on for various 



aspect ratios flxed by the diameters of both cylinders |TU| , pD|^p2| . Also, many experiments 
and comparisons with the theoretical predictions were done |pll- p5| , ^ , p9| , ^lH35| , ^ , ^ | . 
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Our paper is organized as follows: In Sec. |I| we describe the system, the subregions 
of absolute and convective instability, and our methods of investigation. Furthermore, we 
recapitulate the GLE. The next section presents the spatiotemporal behavior of PV patterns 
obtained numerically for two different boundary conditions. In Sec. ^ we analytically and 
numerically elucidate the pattern selection observed within the GLE and the NSE. We 
compare the results with each other, with front propagation, and with experiments. The 
last section gives a conclusion. 



II. THE SYSTEM 

We investigate time- dependent, rotationally symmetric vortex structures in a Taylor- 
Couette apparatus with an externally enforced axial flow. The viscous, incompressible fluid 
is confined to the annulus between two concentric cylinders of inner radius ri and outer 
radius r2. The setup is characterized by two geometric parameters: the radius ratio V = 
and the aspect ratio F, i.e., the quotient of the axial extension of the annulus and the gap 
width d = r2 — ri. Mostly we have used in our numerical simulations an aspect ratio of 
F = 50 and a radius ratio of ?7 = 0.75. The outer cylinder is always kept at rest, while 
the inner one has a rotation rate Qcyi- In addition we impose a small through-flow in axial 
direction. The boundary conditions at ri and r2 were always no slip. The conditions at the 
two ends z = and z = F of the annulus are explained in Sec. |ITT| . 

The flow pattern is described by the momentum balance equation for the velocity field 
u, the Navier-Stokes equations 

((9t + u- V)u= --Vp+z/V^u, (2.1a) 

Q 

and the continuity equation 

V-u = (2.1b) 
which reflects the incompressibility of the fluid. Here u is the kinematic viscosity, g the 
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mass density, and p the pressure. The system is characterized by two dimensionless control 
parameters. The Taylor number 

T=^^ (2.2) 

1 — 1] 

is given by the squared rotation rate Qcyi of the inner cylinder. The Reynolds number 

Re = ^ (2.3) 

V 

is proportional to the mean axial through-flow velocity w. For an axially uniform system 



the homogeneous basic flow state |45 



U (r) = VccF (r) + Wapf (r) (2.4) 



is a linear superposition of circular Couette flow (CCF) 



VccF (r) = + - (2.5) 
r 

in azimuthal direction and of annular Poiseuille flow (APF) 

, , + C In r + D , , 

Wapf (r) = Re (2.6) 

in axial direction. We scale lenghts by c?, times by the radial diffusion time azimuthal 
velocities by the velocity of the inner cylinder fic^^ri, radial and axial velocities by ^. Then 

A = ^ , B = ^ (2.7a) 

1 + n I 1 

C = ^—^-—, D = Cln(l~r]) (2.7b) 

l-r/lnr^' ^ (i _ ^ ^ 



2 



(2.7c) 



Here, 4Re/E = d^p is the dimensionless axial pressure gradient driving the APF. 

At the critical Taylor number Tc (Re) p7| , p!9| , |20| , |31| , which depends on the through-flow 
rate, the basic flow becomes unstable to rotationally symmetric, axially extended PV per- 
turbations via an oscillatory instability. There, a nonlinear PV solution branches off the 
basic flow in an axially inflnite system. We consider the deviation 
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u (r, z; t) = uBr + ve^ + we^ (2.8) 



of the velocity field from the basic flow ( |2.4| ) as the order-parameter field to characterize the 
secondary PV structure. We use the relative control parameter 

corresponding to 

to measure the distance from the onset of PV flow for Re ^ and of stationary Taylor 
vortex flow for i?e = 0, respectively. In this notation 

/i, = and Q(i?e)= J^-j^"^ -1 (2.10) 

Tc [Re = 0) 

is the critical threshold for onset of PV flow. The relation between /i and e is 

f^= , ^ ' V (2-11) 
1 + ec [Re) 

The shear forces associated with the axial through-flow sligthly stabilize the homogeneous 



basic state, so ec (-Re) slightly increases with Re [p^,pT|-|33| . For similar reasons a lateral 



Poiseuille shear flow suppresses the onset of convection rolls with axes perpendicular to the 
flow in a Benard setup of a fluid layer heated from below [^ , |T^] . 

A. Ginzburg-Landau description 

Close to the bifurcation threshold Tc (Re) of PV flow, i.e., for small /x, the flow has the 
form of a harmonic wave, e.g., 

w (r, z;t) = A {z, t) e^^^-^— (r) + c.c. (2.12) 

with a complex amplitude A{z, t) that is slowly varying in z and t. The critical wave number 
fcc, frequency Uc |1T9| , ^ , ^ , |3TH 33[| , and eigenfunction w{r) [|l^,|3] appearing in Eq. (|2.12|) have 
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been obtained from a linear stability analysis of the basic flow state as functions of Re. The 
complex vortex amplitude A {z, t) is given by the solution of the ID complex GLE 

To (a + VgA') = /i (1 + 2Co) A + ^l{l + ici) A" 

-^{l + ic2)\A\^ A. (2.13) 

Dot and primes denote temporal and spatial derivatives in the z-coordinate, respectively. All 
coefficients of the GLE have been calculated as functions of Re for several radius ratios 
7]. As a consequence of the system's invariance under the combined symmetry operation 
{z —z, Re —Re} the coefficients tq, 7 are even in Re while the group velocity Vg 



and the imaginary parts Cq, Ci, C2 are odd in Re |ll5| , |27 | 



It should be noted that the control parameter range of /i over which (|2.12|) gives an 



accurate description of the full velocity field of PV flow, say, on a percent level is indeed 
very small: The asymmetry between radial in- and outflow intensities rapidly grows with 
/i and causes higher axial Fourier contributions ~ e^^^^ |]48|-|50| to the velocity fleld that 
are discarded in the yU — > asymptotics of the GLE approximation ( |2.12| ). However, the 
modulus of the flrst Fourier mode of the vortex structures agrees for Re = as well as for 
Re 7^ quite well with the one predicted by the GLE - cf. Sec. On the other hand, the 
PV structure selected according to the GLE differs from the one resulting from the full fleld 



equations - cf. Sec. |V. 



B. Absolute and convective instability 

For small e and Re the control parameter plane is divided into three stability regimes - 
cf. Fig. |l| - characterized by different growth behavior of linear perturbations of the basic 
flow state. Below the critical threshold ec (-Re) for onset of PV flow (dashed line in Fig. |I]) 
any perturbation, spatially localized as well as extended, decays. This is the parameter 
regime of absolute stability of the basic state. 

Perturbations of the basic state can grow only for e > edRe). However, in the presence 



of through-flow one has to distinguish |I3 between the spatiotemporal growth behavior of 
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spatially localized perturbations and of spatially extended ones. The latter having a form 
~ e*'^^ can grow above ec (Re) - in fact ec is determined as the stability boundary of the basic 
state against extended harmonic perturbations. On the other hand, a spatially localized 
perturbation, i.e., a wave packet of plane wave perturbations is advected in the so-called 
convectively unstable parameter regime faster downstream than it grows - while growing in 



the comoving frame it moves out of the system |17| , p!8| , |36| , |51| . Thus, the downstream as well 
as the upstream facing intensity front of the vortex packet move into the same direction, 
namely, downstream. In this regime the flow pattern that results from a spatially localized 
source, which generates perturbations for a limited time only, is blown out of any system of 
finite length and the basic state is reestablished. 

In the absolutely unstable regime (shaded region in Fig. [l|) a localized perturbation 
grows not only in downstream direction but it grows and spatially expands also in upstream 
direction until the upstream propagating front encounters the inlet in a finite system. The 
final pattern resulting in such a situation shows in downstream direction a characteristic 
axial intensity profile under which the PV fiow develops with increasing distance from the 
inlet. 

Within the framework of the amplitude equation the boundary (full line in Fig. |^) be- 
tween absolute and convective instability is given by 



corresponding to e^o„^ = ec + (1 + €c)ficonv Thus the absolutely unstable regime (shaded 
region in Fig. P is characterized by ;U > fi'^^nv o'^; equivalently, by the reduced group velocity 



Vg = ^° Vg = 2J^^ (2.15) 



being smaller than 2. 

It should be mentioned that the GLE approximation ( |2.14| ) of ficonv describes the bound 



ary between absolute and convective instability resulting from the NSE [ P5| , PU | very well for 
the small Reynolds numbers considered here. 
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C. Methods of investigation 

The linear growth analysis of modes exp [i{kz — oot + rrup)] shows that for small through- 
flow rates the homogeneous basic state becomes first unstable to axisymmetric PV flow 



|T9| , |20|j5^ . Up to -Re = 4 these patterns are also detected experimentally [p5|, p8| , p^ ,pl|-|33 
However, at larger Re one observes stationary spirals and mixed flow patterns [^,^ in 
addition to the bifurcation of propagating spirals P3|,PT]-p3H^ . The latter are predicted 



by the linear stability analysis to branch off the basic state at Re ~ 20 |T^,|2D|,|^ . The 



PV patterns occuring at small through-flow rates {Re < 4) that are discussed here are 
rotationally symmetric. Therefore it is sufficient to solve the hydrodynamic field equations 
in an r-z cross section of the annulus to describe the resulting field u(r, z; t). 

We have performed numerical simulations of the 2D NSE. They are compared with an- 
alytical and numerical results obtained from the ID GLE. The latter was solved with a 
Cranck-Nicholson algorithm using central differences for spatial derivatives with a resolu- 
tion of 20 grid points per unit length d. The solution of the NSE was obtained with a 
time-dependent finite-differences marker and cell (MAC) algorithm |]6|j53| with pressure and 



velocity being iteratively adapted to each other with the method of artificial compressibility 
Also here the spatial resolution was 20 grid points per unit length d. The temporal 



stepsize was 1/1800 times the radial diffusion time d'^/u. 

When comparing finite-differences solutions of the NSE with experiments or with ana- 
lytical properties, e.g., of the GLE we take into account that the critical properties of the 
finite-differences MAC code differ slightly from the latter due to its finite spatiotemporal 
resolution. In particular the critical Taylor number, Tc{Re), of the MAC code lies slightly 



(less than 1.5%) below the theoretical bifurcation threshold - cf. Table I and p^. To find 
the marginal stability curve Tstab{k) of the MAC algorithm for axially extended PV pertur- 
bations of wave number k we analyzed the complex growth rate s{k) in systems of length 
27r//c using periodic boundary conditions in axial direction for through-flow rates up to 
Re = 5. From this analysis we also obtained the critical values of the frequency Uc, group 



February 9, 2008 12 

velocity Vg, wave number kc, and the parameters Tq and (Table I). Furthermore, we have 
cross-checked these results by investigating the evolution of localized perturbations in long 
systems of lengths F > 50, in particular in the convectively unstable regime in a manner 



that is quite similar to experimental procedures p3| , p9| : We generated tiny, localized PV 
perturbations of about five vortex pairs with wave numbers close to kc under an intensity 
envelope of Gaussian shape. They propagate downstream with the group velocity Vg. A fit 
to the rate of exponential growth of the envelope gives /i/ro. From the growth of the full 
width at half maximum we obtain ^q/tq. These identifications are based upon comparing 



with the temporal evolution of a Gaussian perturbation 

v2 



A{z,t = Q)^exp[-^] (2.16) 



of initial axial extension Aq according to the linearized GLE. It yields 

^(^, « ^ exp [/i(l + ico)^^ - ] (2.17) 

with 

^\t) = /\l + 2e,{l + ic,)-. (2.18) 

To 

Test runs with twice the lattice points per unit length showed that the differences between 



the critical properties of the continuous system that were obtained |2^ with a shooting 
method and the finite-differences system significantly decrease by one order of magnitude. 

We have also investigated the dependence of nonlinear, saturated PV flow on the nu- 
merical discretization. These tests show that the nonlinear vortex structures are basically 
independent of the discretization - provided it is not to coarse. However, one should base 
the comparison of bifurcated flow structures obtained with different discretizations on the 
relative control parameter ji = ^ — 1 that is influenced via Tc by the discretization in 
question - see also Sec. |1V B 3 . 



When comparing results obtained for different yU and Re from numerical simulations of 
the NSE with those following from the GLE we found it sometimes advantageous to present 
them as functions of the reduced group velocity Vg (|2.15| ). The deviation 
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2 - l/^ = 2 - 2v//i^_//x (2.19) 

is a scaled distance from the boundary between the absolutely and convectively unstable 
regime. Its relation to the through-flow Reynolds number Re and the Taylor number T = 
(1 + fi)Tc{Re) can be read off from Table I. 

III. PROPAGATING VORTEX PATTERNS 

We investigate PV patterns in the absolutely unstable regime in an annulus of finite 
length with two different end conditions: a basic state boundary condition (BCI) in Sec. [Ill A 



and an Ekman vortex boundary condition (BCII) in Sec. |111 B|. For both boundary condi 



tions numerical simulations have been performed for the parameter combinations marked 
by symbols in Fig. |I|. In each case we observe a pattern of vortices propagating downstream 
under a stationary intensity envelope after transients have died out as shown in Fig. ^. The 
oscillation frequency u of the flow is independent of the radial and axial position while the 
local wave number k, the phase velocity Vp = u/k, and the vortex flow intensity varies 
with r and z. The frequency and the spatial variation of the PV pattern depends on the 
control parameters and on the boundary conditions but not on parameter history or initial 
conditions. 

A. Basic state boundary condition — BCI 

Here we discuss vortex suppressing boundary conditions that are realized by imposing 
the homogeneous basic state U(r) (|2.4| ) at in- and outlet of the annulus. The flow resulting 
for this boundary condition resembles the experimental one of Babcock et al. |33|. There, 



a system of flow distributors and meshes at the inlet reduces external perturbations pen- 
etrating the interior of the annulus and seems also to suppress radial flow Near the 
inlet no stationary Ekman vortices are visible in this experiment and the amplitude of the 



downstream propagating vortices starts with zero or nearly zero ||3T|-p3|. 
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To illustrate the global properties of the flow patterns we present in the upper part 
of Fig. a hidden-line plot of the axial velocity field w. Thin lines show snapshots of w 
at the radial position r = ri + 0.225 obtained from the NSE at successive, equidistantly 
spaced times after transients have died out. Then vortices propagate downstream under 
a stationary intensity envelope (thick line). This envelope is determined by the temporal 
extrema of w{ri + 0.225; t) at any z-position. 

The oscillation frequency of the PV pattern is constant over the entire annulus. Suffi- 
ciently away from inlet and outlet we observe a bulk region of nonlinear saturated PV flow 
with spatially uniform amplitude and wavelength. With increasing through-flow this bulk 
pattern is pushed further and further downstream. The growth length / from the inlet over 
which the amplitude reaches half its saturation value depends on fi and Re. It increases 
and flnally diverges when the control parameters /i, Re approach the absolute-convective 
instability border /i^"™, that is when Vg ( p.l5| ) reaches 2. 

In Fig. ^ we compare the scaled growth length 

L = ^l/io (3.1) 

computed from the amplitude equation (solid line) with results from the NSE (symbols) for 
various combinations of e and Re. Due to the scaling property of the GLE ||T^-[l6| keeping 
in mind the smallness of the imaginary parts q all values for / obtained by the GLE subject 
to BCI fall onto one curve in the plot of L vs. Vg. The open symbols representing the NSE 
results for BCI lie very close to the GLE curve. Full symbols for ECU are discussed in 



Sec. IIIB 



For further characterization of the PV flow structure we found a temporal Fourier de- 
composition of the time-periodic flelds, e.g., 



w r, z. 



t)=Y.Wn{r,z)e-'''^' (3.2) 



to be useful. To that end we flrst determined the oscillation frequency uo of the PV pattern 
and checked that it showed no spatial variation. In Fig. ^ the zeroth temporal Fourier mode 
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Wq and the moduli \wi\, \w2\, and {w^l resulting from the NSE are shown by full lines while 
\wi \ from the GLE is shown by the dashed one. The latter was rescaled by a factor of ~ 0.94 
which comes from comparing the first axial Fourier mode of w{r = ri + 0.225) at Re = in 
an axially periodic system with the GLE result of Recktenwald et al. [0,0 . As an aside we 
mention that for the radial velocity u having a nodeless radial profile with a single maximum 
the deviation of the first Fourier mode of u{r = ri + 0.5) between finite-differences NSE and 
GLE is only 0.8% at Re = 0. 

Note that the solution of the GLE (|1|) 

A(z,t) = /?(z)e'[^(^)-™] (3.3) 

oscillates harmonically with frequency Q under a stationary envelope 

Riz) = \Aiz,t)\ (3.4) 

after transients have died out. Therefore, the GLE velocity field ( |2.12| ) contains no temporal 
Fourier mode other than n = 1 whereas the solution of the full NSE has higher harmonics. 

In the bulk region the temporal modes obtained for finite through-flow from the NSE 
grow for small /i proportional to /i"/^ with relative corrections proportional to /i. Thus, 



they show the same growth behavior with /i that the axial Fourier modes i48|j57|j58|J49|j50|l of 
stationary Taylor vortices without through-flow show as a function of e. The reason is that 
all fields in the PV state have the form of propagating waves 

Mr,z;t) = Mr,z-^t) (3.5) 

in the bulk region - denoted by a subscript b. There the wave number kb and the phase 
velocity uj/kb do not vary with z. 

The Fourier mode Wq of the PV pattern is not zero everywhere since the variation of the 
flow amplitude under the fronts causes a secondary, stationary, closed flow pattern there - 
see, e.g., at 2 ~ 9 and near the outlet in Fig. |a. Its strength and structure depends 
on the steepness of the front. For example, an increase of Re leads to a smaller gradient 
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of the upstream facing front and to weaker stationary secondary flow extending further 
downstream. 

The intensity variation (cf. Figs. 0, ^) of PV flow under the fronts also causes there a 
spatial variation in the local wavelength X{z) (Fig. ^) and in the phase velocity Vp{z) = 
uj\{z)/2'K with the oscillation frequency uj of PV flow being constant. The full (dashed) 
line in Fig. is the wavelength profile \{z) selected for BCI according to the NSE (GLE). 
This figure shows that the GLE does not describe the pattern selection quantitatively. The 
selected wavelength in the bulk region of the NSE (GLE) is smaller (larger) than the critical 
one - see also Sec. |^ 

The local wavelength shown in Fig. ^3 was determined from the phase gradient of the 
first temporal Fourier mode 'Wi{z) at the radial position r = ri + 0.225. We have analyzed Wi 
and also Ui at other r positions and in addition we determined the axial distances between 
node positions of u and w. Under the fronts there is a significant radial variation of A. 
Furthermore, the wavelengths determined via the phase gradients of Wi and Ui differ there 
from each other and from those obtained via node distances of w and u. But in the bulk 
part of the PV structure all these quantities yield the same wavelength A;,. 



B. Ekman vortex boundary condition — BCII 

In Taylor-Couette experiments without through-flow rigid nonrotating end plates bound- 
ing the fluid in axial direction induce stationary Ekman vortex flow [p|,^,|59|. These vortices 
are also detected by Tsameret et al. in experiments with axial through-flow, where meshes 
at the apertures were used as nonrotating boundaries at in- and outlet p8|j29[ . 



The second axial boundary condition BCII enforces stationary Ekman vortices near the 
apertures of the annulus. To that end we impose at both ends 2; = 0, F zero radial flow, 
zero azimuthal flow, and in axial direction the annular Poiseuille flow WAPF{f) ( P-6| ). The 
spatiotemporal properties of the vortex pattern subject to this BCII can be seen in the lower 
part of Fig. There we show a hidden-line plot of the axial velocity in the same way and 
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for the same parameters as for the BCI pattern in the upper part of Fig. |[ 

In the immediate vicinity of in- and outlet there are stationary Ekman vortices whose 
intensity rapidly decreases towards the bulk of the annulus. In addition there is - as in the 
BCI case - the PV flow structure that oscillates in time with a z- independent frequency lj. 
The oscillation amplitude drops to zero near the apertures due to the boundary conditions, 
and the stationary Ekman vortices increase the growth length / of the oscillating structure 
in comparison to the BCI case. 

These two different flow elements are best separated by a temporal Fourier analysis. The 
results for BCII are shown in Fig. |^ for different through-flow rates. The zeroth temporal 
mode Wo (thick line) represents the stationary Ekman vortex flow of the system. The 
temporally oscillating PV structure is characterized by the modes w„ (thin lines). The 
Ekman vortex structure at the inlet is only slightly affected by the through-flow: the strength 
of the Ekman vortex closest to the inlet decreases somewhat with increasing Re, while its 
extension slightly increases. The stationary vortices at the outlet become more and more 
squeezed together with increasing Re, and also their intensity reduces. On the other hand, 
the oscillatory pattern of PV responds dramatically to the through-flow in being more and 
more pushed downstream. Note also that for small Re where the PV amplitude of, say, \wi\ 
overlaps near the inlet with the Ekman vortex intensity the latter causes axial oscillations 
in the harmonics of the PV flow and similarly near the outlet - cf. Fig. 

The filled symbols in Fig. ^ represent the growth length of the PV flow intensity. This 
length diverges in a similar way at the convective instability border = 2 as the one 
obtained for BCI. But due to the presence of Ekman vortices at the inlet the characteristic 
growth length /(e. Re) of PV flow is increased by 0.5 — 2.5d in comparison to the BCI case 
- Ekman vortices push the PV structure further downstream. The difference Ibcii ~ Ibci 
decreases with increasing through-flow rate as Ekman vortices and PV flow become more 
and more separated from each other. 

As in the BCI case the zeroth temporal Fourier mode exhibits a secondary, stationary, 
closed flow pattern of low intensity under the fronts of the PV pattern. For small Reynolds 



February 9, 2008 18 

numbers this pattern is concealed by the Ekman vortex flow, whereas it becomes visible for 
bigger through-flow rates. 



IV. PATTERN SELECTION 

Our investigations of the PV structures show that the axial through-flow causes a unique 
pattern selection such that the selected structure is independent of history and initial con- 
ditions and depends only on the control parameters and boundary conditions. Thus, the 
axial flow causes possible wave numbers within the Eckhaus-stable band of stationary Taylor 
vortex patterns without through-flow to collapse to only one uniquely selected PV structure 
for nonvanishing through-flow rates. However, the PV patterns that are selected according 
to the GLE approximation and to the full NSE differ distinctly from each other. We flrst 
investigate in Sec. [IV A| the selection mechanism within the GLE framework. In Sec. [IV B 



we compare with NSE results and experiments. 

A. Pattern selection within the GLE 

Here we elucidate how the pattern selection mechanism of the GLE can be understood 
as a nonlinear eigenvalue/boundary- value problem where the frequency of the PV pattern 
is the eigenvalue. Thus, the selected PV structure is characterized and determined by the 
combination of eigenvalue and corresponding eigenfunction. 

1. The eigenvalue problem 
We look for solutions of the GLE ([2l3| ) of the form (U) 



A{z, t) = a{z)e-'^' = R{z)e'^^'^'^-^'^ (4.1a) 
with stationary envelope R{z), stationary wave number 

q{z)=k{z)-K = ip'{z), (4.1b) 
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and constant frequency 

Q = u — ujc- (4-lc) 



Inserting this solution ansatz into the GLE (|2.13|) one obtains the nonhnear eigenvalue 
problem 

To {—iila + Vgu') = yU (1 + ico) a + (1 + ici) a" 

--f {1 + ic2) \af a (4.1d) 



or, equivalently, 



iro(-f2 + Vgq)R + T^VgU! = + ico)R 
+eo(l + ic{){R!' - q^R + iq'R + 2iqR') - 7(1 + ica)/?^ 



(4.1e) 



as a solvability condition with Vt being the eigenvalue. We are interested in solutions R{z) 
and q{z) that look like the dashed lines in Fig. i.e., where the variation of R{z) and in 
particular of q{z) is "as small as possible". Such a solution type seems to be connected 
with the eigenvalue VL that is closest to zero. Here we consider the boundary conditions 
Rin,out = at in- and outlet. These requirements fix four boundary conditions - namely, 
Re{Ain,out) = Ifn{Ain,out) = - which are necessary to solve the GLE (|2.13| ). 

Note that the eigenfrequency Q is in general determined by global properties and not 
by the local variation of the eigenfunctions R{z) and q{z), say, at the inlet or in the bulk. 
Nevertheless it is informative and useful for our further discussion to list below relations 
between Q and structural properties at the inlet and in the bulk of the annulus. 

2. Inlet behavior 



From ( [4.1e| ) one finds that the BCI at in- and outlet, R = 0, implies the relation 



Rin,out ~ if2/i I \ '^^'iin,out\Rin,out (4-2) 



February 9, 2008 20 

at z = 0, r. Hence, whenever the modulus grows with finite slope, 7^ 0, the boundary 

condition R = fixes the wave number at inlet and outlet to the value = Qout 

c _ 1 '^^'"9 _ Cl 

^ 4oli + Cij \/l + Ci 



that depends only on Re. This value ensures that the imaginary part of the square bracket 
in ( |4.2| ) is zero. For _Re = one obtains qin,out = as Cross et al. [Q. In addition the 
condition i? = yields via the real part of (|4.2|) to two relations 



between the different slopes 7^ -R^^j and curvatures 7^ -R^'^^ of the modulus at the 
two respective boundaries that hold independently of /i. 

Furthermore, one can derive a relation between the eigenfrequency f2, the wave number 
(lin,out, and its slope q'inout i^l^t and outlet. To that end we consider the axial derivative 
of ( [4.1e|) at inlet and outlet with R = 

iTo{n - VgO) + yu(l + icq) - ^0(1 + ici)q^ 
R" R'" 

(4.5) 

This allows to solve for and q' separately. Using ( ^.3| , |4.4| ) one obtains the relation 



TqVL = (ci - Co)/i + roVgqin,out 

-m + cl)eoqL,ouf (4.6) 

Note that this equation implies q^^ = q'^^^. 

The relation ( [4.6|) also demonstrates the nonlocality of the pattern selection mechanism 
and of the eigenvalue problem. At first sight one might be tempted to infer from the above 
relation between Q and that the latter two fix the frequency. However, Eq. ( |4.6|) 

equally well holds for qout, q'out outlet properties do not fix Vt neither. So the correct 

interpretation of Eq. ( |4.6| ) is that the eigenvalue VL being the characteristic global signature 
of the pattern fixes the local quantities q'ln.out given boundary conditions. 
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3. Bulk behavior 

The wave number q{z) away from in- and outlet in general differs from qin{= Qout)- In 
order to obtain the full axial profiles of the eigenfunctions R{z) and q{z) belonging to the 
eigenvalue fl one has to solve the eigenvalue problem ([4.1|) with the boundary condition 
Rin,out = 0. Alternatively one can solve numerically the time-dependent GLE ( p. 13 ). 



However, when the system size and the control parameters are such that the PV pattern 
forms a bulk part with a homogeneous modulus Rb and wave number qj,, i.e., where 

Ri = R'l = qi = 0. (4.7a) 

Then, the dispersion relation 

Tofi = (C2 - Co)n + ToVgqb + (ci - C2)^Qql (4.7b) 

provides a relation between the frequency eigenvalue fl, the bulk wave number qh, and the 
bulk modulus 

Rl = /iniM. (4.7c) 

7 



The equation ( |4.7b[ ) establishes together with ( |4.6| ) also a relation between qi, on the one 



hand and qin = qout and = g^ut on the other hand. 



4- Pattern selection at the convective instability boundary 



Approaching the convective instability boundary ( p.l4|) we found that q^^ decreases to 
zero so that the frequency eigenvalue becomes 

'^O^'conv = -(CO + Cl)/ieo™. (4.8) 

Then the bulk wave number approaches according to ([4.8| , |4.7b| ) the limiting value 



UqbYcon. = -VJ^., (4.9) 

Cl - C2 
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while the inlet wave number gj„ is given by ( |4.3| ). Due to the very small imaginary parts Cj 
the solution ( [4 .91) with the plus sign has to be discarded, since it corresponds to unphysically 
large wave numbers. 

In Fig. ^ the bulk wavelengths selected according to the GLE are shown by halftone 
symbols. The results in the absolutely unstable regime, Vg < 2, were obtained numerically 
by integrating the time- dependent GLE (|2.13|) . The value shown right at the border, Vg = 2, 
of the absolutely unstable regime is the analytical expression (^]9|). Note that the selected 
GLE wavelengths are mostly larger than the critical one. 



5. Comparison with front propagation 

It is suggestive to compare properties of PV patterns like the one in the upper part of 
Fig. 1^ with those behind an upstream facing front that connects in an axially infinite system 
the basic state at z = — oo with the developed PV state at ^ = +oo. A linear growth 
analysis of the GLE ( p.l3| ) along the lines of Sec. VI B 3] shows that the far tail of a 
linear front of PV perturbations moves with the front velocity 

VF = Vg-2^Jfi{l + cl). (4.10) 

In the convectively unstable regime < fi < ^i^onv ^^e front propagates downstream so that 
vp > 0. In the absolutely unstable regime fi > fJ^conv moves upstream, i.e., vp < 0. And 
at the boundary /i = /i^o^^ the front is stationary - vp = 0. We compare the properties of 
such a stationary front at /i = fi^onv "with those of the solution A{z,t) ( |3.3| ) with stationary 
modulus in a semi-infinite system (0 < 2; < 00). The reason for restricting the comparison 
to jJ^conv is that an upstream moving front in the absolutely unstable regime will be pushed 
against the inlet whence the spatiotemporal structure of the free front gets modified by the 
inlet. 

The far tail of the stationary front resulting from the linear GLE is characterized by the 
local wave number 
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^oqpivF = 0) = j^L=./J]^ (4.11) 

1 + cf 



and the oscillation frequency 

TonpivF = 0) = -(co + Ci)lJ,l^^^. (4.12) 
In fact in the linear part of the stationary front the flow amplitude 

A(2,t) ~e'['3»"-C(Q»)tl (4.13) 

varies with a wave number Re{Qs) = qp ( |4.11|) , a spatial decay rate Im{Qs), a temporal 

growth rate Im{Q{Qs)) = 0, and an oscillation frequency Re{({Qs)) = ( [1-12|) . Here 

^oQs = -Iri-rr^, = -4i=v& (4.14) 



is the saddle position of the complex dispersion relation 

roC(<5, = T-ot^sQ + ^(1 + «co)/i - z(l + ici)^oQ^ (4-15) 



of the linear GLE in the complex Q-plane determined by the condition 

dCiQ) 



dQ 



= 

Qs 



(4.16) 



for vf = 0, i.e., at fi = 

So the inlet wave number gj„ ( [4.3|) agrees with qp ( [4.11j ) and ficon?; coincides with 
(|4.12|) . Thus, right at the border fi'^onv '^^ the absolutely unstable regime the frequency 
eigenvalue of the PV pattern growing in downstream direction from the inlet at 2; = agrees 
with the oscillation frequency selected by the stationary, linear front. Within the front 
propagation point of view the PV pattern develops from the basic state with an exponential 
growth rate k, of the modulus. This leads for R ^ to R' = kR,R" = k^R, ... — > for 
z ^ —00. On the other hand, the PV flow intensity in the semi- infinite system in general 
drops to zero at the inlet with finite R' and R". Only at the border line /i = ^^onv 
difference disappears since there R', R", ... ^ for ^ = 0. In the absolutely unstable regime, 
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Vg < 2, the frequency eigenvalues Q = u — Uc of the GLE he shghtly above the hmiting 
value flconv the border of the absolutely unstable regime. The shaded symbols in Fig. |^ 
show for three different e how Q obtained numerically by integrating Eq. ( p. 131 ) approaches 
with increasing Re the front frequency Re{({Qs)) = ^f{vf = 0) = ^^onv (dashed line). 
For convenience, the merging points with this limit line are marked for the three e-values 
investigated here by small halftone symbols. 

Using this frequency for the full nonlinear PV pattern at f^conv obtains from the bulk 
dispersion relation ( [4.7b|) an expression for the bulk wave number 



c 

conv 



VJ^v (4.17) 



Cl - C2 

far behind the front. 

This bulk wave number behind the front has the opposite sign of the one that is cited 
in [^] and Eqs. (10,14)]. An explanation for this might be that the formulas of 



Nozaki and Bekki [^, Eqs. 2, 8, and 9] that have been used in refs. [^,^ do not seem to 



have been transformed to a through-flow situation with an upstream facing front that joins 
in upstream direction to the basic state, in our case at z = — oo. After correcting the sign 
error, the rough qualitative agreement that was reported in [^,^] between the experimental 
wave numbers [0,^ and the erroneous GLE results [p8| , p9| disappears. In fact, both, the 



experiments [28,29| and our numerical simulations of the NSE yield wave numbers that differ 



in a common distinctive way from the proper GLE result ( |4.17|) . 



As an aside we mention, that the nonlinear front solution of Nozaki and Bekki pi 



Eqs. 2-7] yields for our GLE ( p.l3| ) a front that is not stationary at ficonv but rather moves 
with velocity 



vnb = -(3, 



^ 8 + 9c? 

in upstream direction with a bulk wave number 



^ + '' IK (4-18) 
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b NB 



{c2-ci)J8 + 9ci 



X 



3(1 + cf) + sign{c2 - Ci)J8{c2 - Ci)2 + 9(1 + ciCs 



(4.19) 



So we conclude that the nonlinear front solution of Nozaki and Bekki is unrelated to the 
PV patterns at /i^ 



^conv ■ 



Furthermore, it was noted |£2| that localized initial perturbations did 
not evolve into the nonlinear front solution of Nozaki and Bekki |^.B^,B3|. However, vnb 
and {qbjNB differ for small through-flow only slightly from the respective values Vp = ^ and 
[qb)F- For example at Re = 2 one obtains vnb = —O.OGvg and {qbjNB = l-15(g6)F- 



6. Scaling properties 



Scaling length and times according to 



z = ^^ ; i = fi-, (4.20) 
^0 To 



the GLE for the reduced amplitude A = A^^/ ^ does no longer contain 7 and ^ explicitly 
but the reduced group velocity Vg ( p.l5| ) and the coefficients cq, ci,C2. Thus the reduced 



selected frequency VL = tqVL/^ and the selected bulk wave number = ^oQ'b/v^ depends 
not only on Vg but via Co,Ci,C2 also on the Reynolds number. This -Re-dependence can 
alternatively be seen - via the i?e-dependence of Vg and of the q's entering ( |2.15| ) - also 



as an additional /i-dependence. The latter dependence is sufficiently strong to prevent a 
scaling of f2 with Vg alone. However, by reducing the eigenfrequencies f2(e. Re) selected for 
our three different e-values by the limiting values fi(e, Re^^^^^) = ^conv (N^-8|) at the border 
of the absolutely unstable regime we effectively have eliminated the e-dependence and all 
GLE data (halftone symbols in Fig. |^) fall onto one curve. A similar scaling holds for the 
bulk wave number divided by {qb)conv (f^-QP- 

Remarkably enough, by reducing the NSE frequencies with the GLE frequency ^conv 
( [4.8| ) the NSE results (open symbols in Fig. P) almost show this one-variable scaling with 
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Vg. However, for the GLE (cu — c<Jc)/^^coni; increases monotonously with the scaled group 
velocity Vg whereas the NSE results show a monotonous decrease with Vg. This discrepancy 
between NSE and GLE results is caused by the different dispersion relations of the full 
hydrodynamic equations and the approximate GLE - cf. below. 

B. Pattern selection within the NSE 

Here we compare our numerical solutions of the NSE with results from the GLE and 
from experiments. First we discuss the selected wave numbers in the bulk region of PV flow. 
Then in comparison with front propagation we show common selection properties of NSE 
and GLE. Finally we compare with experiments. 

1. Bulk wavelengths - NSE vs GLE 

In Fig. ^ we show all selected bulk wavelengths obtained from numerical simulations 
of the NSE for BCI (open symbols) and BCII (filled symbols) in comparison with results 
obtained numerically from the GLE for BCI (halftone symbols) . The NSE wavelengths that 
are selected when the basic state is enforced at in- and outlet weakly decrease with increasing 
Re and e. For the smallest e = 0.0288 their values are very close to the critical wavelength 
Ac (dotted line). For this boundary condition the wavelengths are mostly smaller than the 
critical one with a deviation of up to 2% for the largest e = 0.186 that was investigated 
here. These deviations are much stronger than the decrease of the critical wavelength with 
increasing Reynolds numbers. 

It is interesting to note that transversal Rayleigh-Benard convection rolls in horizontal 



shear flow p!^ , p!5| show a qualitative and quantitative [p^ e = 0.114] correspondence to 
PV flow. Therefore, our results are not limited to PV flow in a Taylor-Couette apparatus 
only. In this context we also mention the experimental investigation of Gu and Fahidy 
PB| . They observed that axial through-flow causes a rearrangement of the vortex centers 
similar to the one seen for tranverse convection rolls subject to a horizontal flow |]15|, Fig. 2]. 
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Without through-flow the vortices are ahgned along a straight hue roughly in the gap center. 
With increasing through-flow vortices are alternatingly displaced towards the outer or inner 
cylinder, thereby reducing the axial component of their distance and thus the wavelength. 
Viewed in the z — r plane with r2 above ri, right (left) turning vortices - for which w 
is negative (positive) close to the inner cylinder - are displaced towards the inner (outer) 
cylinder by a through-flow directed from left to right, since the thro ugh- flow enhances the 
axial flow component w of a right turning vortex near the outer cylinder at r2, and weakens 
it near the inner one at ri and vice versa for the left turning one. 

Obviously the wavelengths resulting for BCI from the NSE (open symbols in Fig. ^) differ 
substantially and systematically from those resulting from the GLE (halftone symbols): the 
former decrease with growing Re and e while the latter increase. Nevertheless we think - 
cf. the discussion in Sec. |1VB 4| - that the selection mechanism is quite similar. 



2. The effect of Ekman vortices 

The NSE wavelengths obtained for the Ekman vortex generating BCII (flUed symbols 
in Fig. P) are similar to the BCI results in the limiting cases of small Re or small e, and 
near the border line = 2 of the absolutely unstable regime. In the vicinity of = 2 
the stationary Ekman vortices and the PV flow are spatially separated, and their interac- 
tion is very weak allowing for PV wavelengths like those for BCI. For very small Reynolds 
numbers the agreement in the PV wavelengths for BCI and BCII is not yet understood. In 
between there is a visible difference of up to 5% (e = 0.186 and Vg ~ 0.4) between the bulk 
wavelengths obtained for the two boundary conditions. 

However, we found that the bulk phase velocities, Vp = uj/k^ were almost independent of 
the boundary conditions: the relative difference of Vp for BCI and BCII is less than 0.1% for 
the parameter regimes investigated here. Furthermore, the deviation of Vp from the critical 
value (fp)c = ^c/kc is less than about 1%. Hence, one can infer from the plots in Fig. |^ for 
A/ Ac — 1 also the variation of the frequency eigenvalue with through-flow since 
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k 

- 1 ~ — - 1. 



(4.21) 



For the largest e = 0.186 we have observed for the Ekman vortex generating BCII 
an interesting behavior for very small Re < 0.5 (Vg ^ 0.3). For these parameters the 
PV pattern did not develop an axially homogeneous wavelength in the spatial region where 
the flow amplitude was indeed spatially uniform. Under these circumstances the wavelengths 
roughly show a linear variation with z between a larger value (upper small filled triangles 
in Fig. P) at the upstream end of the bulk amplitude region and a smaller value (lower 
small filled triangles in Fig. |^) at the downstream end of the bulk region of constant fiow 
amplitude. So the Ekman vortex at the inlet inpedes the free downstream motion of the PV 
fiow - it exerts a phase pinning force on the PV pattern that stretches the node distances. 
The situation is somewhat similar to the explicitly phase pinning boundary conditions that 
have been investigated in Rayleigh-Benard convection with through- fiow p5|jr5| . 



The selected wavelengths in the bulk region are practically independent of the system 
length as long as a saturation of the Taylor vortices takes place. To check this we performed 
a few simulations with aspect ratios F = 25 (e = 0.114, i?e = 2, BCI and BCII) and 
F = 10 (e = 0.114, 0.186, -Re = 1, BCI and BCII) and compared with the results obtained 
for our standard length F = 50. The change in the selected wavelength or frequency is 
less than 0.1% if the aspect ratio is halved to F = 25. In the system with length F = 10 
no bulk region can be observed, therefore the wavelength can not reach a constant value. 
The oscillation frequency of the fiow pattern is always independent of the axial position 
and mainly unchanged; the deviations are less than 0.5% even for the very short system of 



We have also investigated the dependence of nonlinear, saturated PV fiow on the nu- 
merical discretization. These tests show that the nonlinear vortex structures are basically 
independent of the discretization - provided it is not to coarse - if one bases the comparison 



3. Finite size effects? 



V 



10. 
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of bifurcated flow structures obtained with different discretizations on the relative control 



parameter /x 



Tc 



— 1 that is influenced via by the discretization in question [36 



4- Comparison with front propagation - NSE vs GLE 



Let us first consider Taylor vortices without through-flow. Then the GLE ( |2.13| ) contains 



no imaginary coefficients Cj and the first-order spatial derivative is absent since Vg = 0. In 
this case the GLE yields for propagating fronts ||6^ as well as for stationary patterns in 



finite systems the critical wave number, q{z) = 0, all over the extension of the pattern 
- the boundary condition A = enforces the collaps of the supercritical band of stable 
bulk wave numbers of nonlinear vortex patterns to g = 0. This nonlinearly selected wave 
number q = happens to be for Re = the same as the wave number of maximal linear 
growth under a linear front whose spatiotemporal evolution is governed by the dispersion 
Cgle{Q, e; -Re = 0) ([4.7b| ) of the linear GLE. In the presence of through-flow, however, the 
wave number of maximal growth under the linear front - e.g. qp ( [4.11| ) - of the GLE differs 



from the one in the nonlinear hnlk far behind the front [{qb)F ( |4.17| )]. 

Now, already without through-flow the dispersion relation of the linear NSE differs sig- 
nificantly from that of the GLE so that the wave number of largest temporal growth under 
a linear NSE front deviates from the GLE result g = 0. As an aside we mention that a 
somewhat similar behavior has also been observed experimentally P,p!0| and numerically for 



Taylor vortices or Rayleigh-Benard convection rolls |jTT| which, however, were located 

in the nonlinear part of propagating fronts. In this context it is useful to keep in mind 
that the wavelength profile X{z) of numerically obtained Taylor vortex fronts []6|,|68| shows 
an axial variation with a characteristic dip in the steepest region of the front similar to the 
one of the full line in Fig. ^ near z = 8. 

Now let us consider the through-flow case. In Fig. |^ we compare the eigenfrequencies of 
the nonlinear NSE (open symbols) for BCI with those (solid line) that would be selected 
by a linear front that is stationary at the border of the absolutely unstable regime. To that 
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end Recktenwald and Dressier have determined the dispersion relation ({Q) of linear 
perturbations ~ g«[(fcc+Q)2-('^c+C)t] jj^ the complex Q-plane resulting from the NSE with a 
shooting method as described in ref. |^3[. The solid line in Fig. |^ represents the frequency 
Qp = R^iCiQs)) selected by the stationary front {vp = 0) of the linear NSE. Again, as 
discussed for the GLE in Sec. |1V A| 5, this frequency is determined by the saddle Qs of C(Q) 
for which the temporal growth rate Im ({Qs) = 0. This defines the boundary between 
absolute and convective instability. With increasing Re the eigenfrequency of the nonlinear 
NSE approaches the frequency selected by the linear stationary front at the border of the 
absolutely unstable regime as indicated for the three e-values shown in Fig. ^ The front 
frequencies for these e-values are marked by small open symbols on the curve Re{({Qs)) to 
indicate that the eigenfrequencies of the NSE in the absolutely unstable regime (large open 
symbols in Fig. |^) do indeed end at the right positions on this curve. 

This NSE behavior is very similar to the GLE behavior (halftone symbols and dashed 
line). Note, however, that the linear dispersion relations of NSE and GLE are different: the 
front frequency Re{(NSE{Qs)) — 0.0072_Re + 0.0056-Re^ (solid line in Fig. |^ is positive and 
increases with through-flow while Re{C,GLE{Qs)) — — 0.0002i?e — 0.0062_Re^ (dashed line) is 
negative and decreases with through-flow. This difference in the dispersion relations seems 
to be the major cause for the differences in the selected patterns. 

5. Comparison with experiments 



In early experiments by Snyder or Takeuchi and Jankowski the distinction 



between absolute and convective instability was not yet established. PV patterns were 



observed with wavelengths below the critical one for Reynolds numbers Re < 10 pi| , p2 
However, a pattern selection was not reported. The wavelengths were seen to vary up to 7% 
from run to run for the same control parameter combinations and in some cases 5% over 
the spatial extension. Furthermore, the vortex spacings did not change after the occurrence 



of PV flow when increasing the relative Taylor number up to e = 0.15 P^. 
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Takeuchi and Jankowski obtained wavelengths below the critical one for example 
with a deviation of ~ —2.5% for Re = 9.4, e ~ 0.06, rj = 0.5 that is in the same order of 
magnitude of our NSE results in the absolute instability regime. 



The experimental results of Tsameret and Steinberg p8| , p9| exhibit a unique selection. 



similar to the results presented here. Their wavelengths |^ are always less than the critical 
one for axial boundaries that enforce stationary Ekman vortex flow. However, in contrast to 
our results they observed no e dependence of the selected wavelengths but an Re dependence 
that is much stronger than ours leading to a relative deviation from Ac of up to —10% at Re ~ 



3 for 7] = 0.707 p8|. In addition their phase velocity of PV flow, Vp = 1.055i?e ||2^, is much 
smaller than ours, which is close to the critical one. All wavelengths observed for through- 
flow rates above Re = 1 belong to noise sustained structures in the convectively unstable 



region p^, which leads to deviations from Ac of about 3% when one limits the comparison 



to the regime of absolute instability. Since the e-values at which the measurements of p8|j29 
were performed are not given we have estimated an e- value of ~ 0.025 based on the convective 
instability border at Re = 1. For this e the above mentioned 3% deviation of A from Ac is 
an order of magnitude larger than our NSE results or the results ( [4.9| ) of the GLE. 

On the other hand, the experimental wave number found in for a known, published 
e ~ 0.38 and Re up to 4 can be compared directly with our results to check the predicted 
e and Re dependence. In this case all experimental wavelengths decrease up to 10% 
below the critical one at Re ~ 3.7 whereas we would expect from our simulations and from 



14,15] less than half the deviation from Ac 



V. CONCLUSION 

We have investigated propagating vortex structures in an axial flow for different realistic 
axial boundary conditions in systems of finite length. Within the subregion of absolute 
instabilty a unique pattern selection is observed. The selected PV fiow structures are in- 
dependent of parameter history, initial conditions, and system's length (provide it is large 
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enough to allow for a nonlinear, saturated, homogeneous bulk region). But they depend on 
the axial boundary conditions. For conditions enforcing the basic state at the boundaries 
as well as for conditions that enforce stationary Ekman vortex flow near the boundaries 
the PV pattern is suppressed at the inlet and outlet of the annulus. Then one observes 
a characteristic stationary axial intensity and wave number profile, whereas the oscillation 
frequency of the pattern is constant all over the system. The analysis of the appropriate 
Ginzburg-Landau equation for PV flow shows that the selected frequency of the pattern 
oscillation is the eigenvalue of a nonlinear eigenvalue problem for which the axial variation 
of the corresponding eigenfunction is smooth and as small as possible. The eigenfunction is 
the complex pattern amplitude that characterizes the intensity and wave number profiles of 
the PV structure. The GLE intensity profiles agree well with those of the NSE. However, 
there are characteristic differences in the selected frequencies and wave number profiles of 
the full NSE and the GLE approximation. They are identified to be mainly caused by the 
different dispersion relations of the equations. Approaching the border between absolute 
and convective instability the eigenvalue problem becomes effectively linear and the pattern 
selection mechanism becomes that one of linear front propagation. 
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TABLES 

TABLE I. Bold-faced numbers denote critical values and coefficients that are appropriate for 

the finite-differences code with spatial grid size 0.05 used in our simulations of the NSE and that 
have been used for scaling the NSE results. The methods for determining these numbers are 
described in Sec. II C. The imaginary parts ci, C2, C3 are taken from ref. [27]. 7 depends on the 
normalization of the linear radial eigenfunction u (r) which was chosen to be \u (ri -|- 0.5) | = 15.2. 



a = ao 



1 + ^ {Re/a,y 



To 



7 



ao 2420.23 3.1305 0.144 0.0762 8.06 
02 138.62 252.96 -27.23 -49.68 -12.33 



a = aiRe 1 + -^ {Re/as 



Vn 



l«3| 

Co 



Cl 



C2 



ai 1.20 3.647 1/138 1/40.7 1/287 
as - - -51.8 70.5 -97.8 
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FIGURES 

FIG. 1. Stability domains of the basic flow state ( |2.4| ) in the plane of control parameters. 
Numerical simulations have been performed for the parameters marked by circles (e = 0.0288), 
diamonds (e = 0.114), and triangles (e = 0.186). Dashed line is the critical threshold for onset of 
extended PV flow and full line the boundary ( p.l4| ) between absolute and convective instability. 
The parameters entering ( 2.14| ) - cf. Table I - characterize the finite-differences version of the NSE. 



They have been obtained with the methods described in Sec. II C. 

FIG. 2. Spatiotemporal structure of PV flow. Thin lines show vertically displaced snapshots 
of the axial velocity field w{ri + 0.225, z; t) at successive, equidistantly spaced times. Thick lines 
show the stationary envelopes. The BCI at inlet and outlet suppresses any vortex flow there, while 
the BCII suppresses PV flow but induces stationary Ekman vortices. Parameters are e = 0.114 
and Re = 2.5. 



FIG. 3. Scaled growth length L ( ^ ) of PV structures vs scaled group velocity Vg 



Symbols represent lengths obtained from the NSE for different combinations (cf. Fig. |l|) of Re and 
e = 0.0288 (circles), e = 0.114 (diamonds), and e = 0.186 (triangles). Boundary conditions are BCI 
for open symbols and BCII for filled ones. The line shows the scaling behavior of the GLE subject 
to BCI. 

FIG. 4. Structure of PV flow selected for BCI at inlet and outlet. Shown are the axial variations 
of the temporal Fourier modes (a) and of the wavelength (b) of w{ri + 0.225, z; t). Full lines result 
from numerical simulations of the full NSE. Dashed lines come from the GLE. Parameters are 
e = 0.114 and Re = 2.5. 

FIG. 5. Axial structure of vortex patterns subject to BCII. The zeroth temporal Fourier mode 
vjQ (thick lines) of the axial velocity field w{ri + 0.225, z; t) reflects the stationary Ekman vortices. 
The modes Wn with n > (thin lines) characterize the oscillating PV pattern. Parameters are 
e = 0.114 and Re = 1.0 (a), 2.5 (b), 2.7 (c). 



February 9, 2008 41 

FIG. 6. Selected wavelengths in the bulk region of PV stuctures vs scaled group velocity Vg 
( 2.15| ). Open (halftone) symbols refer to the NSE (GLE) subject to BCI at inlet and outlet. 
Filled symbols refer to the NSE subject to BCII at inlet and outlet. Parameters are e = 0.0288 
(circles), e = 0.114 (diamonds), and e = 0.186 (triangles). The small filled triangles are explained 
in Sec. |VB|2. 

FIG. 7. Selected oscillation frequency of PV flow obtained with BCI from the NSE (open 
symbols) and GLE (halftone symbols) vs Reynolds number. The full (dashed) line is the front 
frequency = ReiCiQs)) of the linearized NSE (GLE) at the border between absolute and 
convective instability. Parameters are e = 0.0288 (circles), e = 0.114 (diamonds), and e = 0.186 
(triangles). The small symbols marking ^Ip for these three e values show that the eigenfrequencies 
of the nonlinear equations in the absolutely unstable regime approach the 0,f limit curve at the 
right places. 

FIG. 8. Frequency shift oj — uJc of PV flow vs scaled group velocity Vg ( glsD for BCI. Both, 
NSE (open symbols) and GLE (halftone symbols) results are scaled by the corresponding GLE 
frequency ^conv (0) the convective instability border for the respective parameters e = 0.0288 
(circles), e = 0.114 (diamonds), and e = 0.186 (triangles). The small symbols have the same 
meaning as explained in Fig. ^. 
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